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Summary . 

The  aim  of  this  project  was  to  undertake  a  statistical 
study  of  Monte-Carlo  methods  in  semi-conductor  device  physics 
using  tools  from  the  modern  theory  of  probability  and 
stochastic  processes.  The  Monte-Carlo  approach  has  been  used 
to  obtain  numerical  estimates  of  parameters  by  simulating  the 
physical  process  on  a  computer  rather  than  attempting  to 
solve  the  integro-dif ferent ial  equations  which  arise.  Thus 
this  method  estimates  averages  and  can  be  applied  to  any 
quantity  expressible  as  such  but  may  require  an  enormous 
number  of  simulations  before  statistically  reliable  results 
are  obtained,  especially  if  the  quantity  if  interest  is 
subject  to  very  low  probability  of  occurrence.  A 
probabilistic  study  of  the  behaviour  of  individual  electrons 
was  undertaken,  complementing  the  Monte-Carlo  approach.  The 
methodology  required  the  formulation  of  some  novel  stochastic 
processes  and  suggested  new  types  of  conditions  for 
ergodicity  of  Markov  processes. 
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Research  Objectives. 

The  objective  of  this  project  was  to  undertake  a 
statistical  study  of  Monte-Carlo  methods  in  semi-conductor 
device  physics  using  tools  from  the  modern  theory  of 
probability  and  stochastic  processes.  The  principle  behind  a 
Monte-Carlo  transport  calculation  consists  of  the  simulation 
on  a  computer  of  the  motion  of  an  electron  as  an  alternating 
sequence  of  free  flights  (drift)  in  an  external  field 
followed  by  a  collision  (scattering)  with  a  phonon  or 
impurity.  Between  collisions  the  electron  obeys,  in  the 
semi-classical  approximation,  Newton's  second  law  h  dk/dt  = 
e  F  where  e  is  the  electron  charge,  h  Planck's  constant,  k 
the  wave  number,  and  F  lL<=  external  field.  The  times  between 
collisions  and  the  final  states  after  collisions  are  random 
quantities  whose  probability  distributions  are  determined  by 
the  electric  field  strength  and  the  transition  functions  of 
the  scattering  events.  These  pairs  of  mechanisms, 
deterministic  drift  and  random  scattering,  are  repeated  and  a 
histogram  is  constructed  counting  the  proportion  of  time  the 
electron  occupies  each  state.  Thus  the  Monte-Carlo  method 
estimates  averages  and  can  be  applied  to  any  quantity 
expressible  as  an  average.  It  is,  in  principle,  capable  of 
giving  accurate  results  but  may  require  an  enormous  number  of 
simulations  before  statistically  reliable  results  are 
obtained,  especially  if  the  quantity  if  interest  involves 
states  having  low  probability  of  occurrence. 

The  stochastic  structure  of  these  simulations  is  quite 
attractive  from  a  statistical  viewpoint,  describing  what  are 
known  as  random  evolutions  and  embracing  a  diversity  of 
topics  from  stochastic  processes.  The  scattering  mechanism, 
for  instance,  is  known  as  a  self-exciting  point  process  in 
the  probability  literature.  u  /as  hoped  that  a 
probabilistic  study  of  the  beh^v’nur  of  individual  electrons 
would  complement  the  Monte-Car_  approach  and  provide  insight 
which  could  not  otherwise  be  obtained.  Conversely  the 
physical  problem  suggested  a  class  of  random  processes  whose 
probabilistic  structure  was  quite  attractive  and  whose 
analysis  would  enrich  the  subject  of  stochastic  processes. 

At  the  start  of  this  work  the  P.I.  intended  to  develop 
close  ties  with  the  Electron  Physics  Laboratory  at  the 
University  of  Michigan  in  order  to  acquire  first  hand  details 
of  the  substantive  issues  and  to  maintain  a  fruitful 
collaboration  ,  in  particular  with  Dr  Peter  Blakey.  However 
Dr  Blakey  left  shortly  after  work  began  and  consequently  the 
research  acquired  a  more  mathematical  flavor  than  originally 
anticipated. 
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The  main  goals  of  this  research  were: 

1.  to  formulate  the  Monte-Carlo  approach  rigorously  as 
stochastic  processes; 

2.  to  prove  limit  theorems  relating  to  steady  state 
behaviour ; 

3.  to  interpret  the  limit  theorems  qualitatively  as 
functions  of  electric  field  strength  and  scattering 
mechanisms  which  could  then  be  interpreted  as 
identifiable  physical  properties  of  semi-conductor 
material ; 

4 .  to  determine  whether  the  speed  of  convergence  in  the 
limit  theorems  could  be  improved  by  biasing  and  if  so 
how  feasible  it  would  be  to  implement  computer  automated 
iasing  in  the  Monte-Carlo  runs. 


Status  of  the  Research 

A  detailed  analysis  has  been  made  of  a  model  motivated 
by  the  momentum  randomizing  acoustic  phonon  scattering  in  the 
(100)  valley  of  Gallium  Arsenide.  This  has  been  extended  to 
processes  having  an  upper  and  a  lower  valley  with  intervalley 
scattering  although  these  results  are  less  complete. 

When  an  electron  in  the  (100)  valley  of  GaAs  undergoes 
acoustic  scattering  the  energy  surface  of  the  final  state  is 
spherical  and  the  momentum  is  uniformly  distributed  on  this 
surface,  unlike,  for  instance,  the  non-uniform  scattering  of 
polar  optical  scattering  in  the  (000)  valley  or  intervalley 
scattering.  The  mathematical  details  of  this  model  are  as 
follows.  There  is  a  stochastic  process  k(t)  representing  the 
momentum  of  a  particle.  The  momentum  increases  linearly  with 
acceleration  a  between  scattering  events.  These  events  occur 
at  random  times  {T^}  determined  by  a  self  exciting  point 

process  according  to  the  infinitesimal  rule 

P [scattering  event  in  time  interval  (t,t+Dt)  |  k(s),  0<s<t  ] 

=  X(k (t) )  Dt  +  o (Dt) . 

The  process  k(t)  is  Markovian  and  its  behaviour  is  completely 
determined  by  the  function  X(k) ,  depending  on  which  k(t)  may 
exhibit  transient,  null  recurrent,  or  ergodic  behaviour. 

While  physical  considerations  demand  that  k(t)  be  ergodic, 
nonetheless  all  three  types  of  behaviour  or  possible 
mathematically.  For  example  if  X  grows  slowly  then  few 
scattering  events  occur  when  k  is  large  and  the  motion  is 
almost  deterministic  to  infinity,  hence  transient. 

To  illustrate  the  nature  of  the  results  obtained,  the 
general  cases  of  which  are  technical  and  are  being  prepared 

for  publication,  let  X(k)  =  constant  =  (1  ^  .  Let  denote  k 


■ 


k  just  prior  to 
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just  after  the  n  scattering  event  and 
the  n^  scattering  event 


Theorem  1 . 


Theorem  2 . 


’heorem  3  . 


V  =  U  +  aT; 
U  =  VZ 


where  equality  and  limits  are  to  be  taken  in  the 
distributional  sense.  In  the  first  equation  of  Theorem  3  U 
and  T  are  independent  on  the  right  side  and  T  is  exponential 
with  mean  m  while  in  the  second  equation  V  and  Z  are 
independent  and  Z  is  uniform  on  (-1,1)  .  The  methodology 
involves  establishing  bounds  of  the  form 

P  [  sup  (n>N)  |P  Z.  |  >  e  ]  <  - 

b(l+b) 


Next  let  S (t )  =  jk(s)  ds . 

0 


Theorem  4 . 


S(t) 

t 


a)l  a .  s 


Theorem  5 . 


S  (t )  ~a[Lt 
(t/H)  1/2 


N  (0,  s  ) 


Theorem  4,  as  ergodic  theorem,  has  a  natural  interpretation 
as  Ohm's  law,  since  the  left  side  represents  the  current,  a 
is  proportional  to  the  electric  field,  and  m  is  a  function 
of  the  medium. 


When  X(k)  =  |k|  analogously  results  hold.  Thus 


a.s.  for  some  constant  K. 


This  displays  non-ohmic  behaviour.  In  general 


g  (a) 


for  some  function  g  which  can  be  described  as  the  solution  of 
a  second  order  ODE. 
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Necessary  and  sufficient  conditions  for  ergodicity  seem 
difficult  to  obtain.  The  usuai  approach  would  be  to  verify  a 
condition  due  to  Doeblin  but  this  condition  is  hard  to  verify 
and  it  may  not  hold.  A  promising  alternative  is  based  on  the 
observation  that  the  transition  densities  p(x,y)  of  many  of 
the  above  models  share  the  following  property: 

there  is  a  real  c  >  0  such  that  if  A  =  {y:y<c}  then 
p(x,A)  >  1/2  uniformly  in  x. 

Intuitively  this  should  cause  the  process  to  return 
ergodically  to  each  state.  The  general  implications  of  such 
a  condition  are  currently  under  investigation. 

Multiple  .cattering  mechanisms  have  also  been  formulated 
by  enlarging  the  state  space  to  a  vector  stochastic  process 
(k(t),  i(t))  where  i  is  an  index  identifying  the  valley. 
Scattering  events  may  be  one  of  two  types  involving  either: 

(a)  change  in  k(t),  no  change  in  i(t); 

(b)  no  change  in  k(t),  change  in  i(t)  (intervalley 
scattering) . 

The  marginal  process  k(t)  is  no  longer  Markov. 

The  results  described  above  are  currently  being 
prepared  for  publication. 


j 


